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Abstract 

We report on a search for the decay K^ -^ n^e'^e' carried out by the 
KTeV/E799 experiment at Fermilab. This decay is expected to have a 
significant CP violating contribution and the measurement of its branch- 
ing ratio could support the CKM mechanism for CP violation or could 
point to new physics. Two events were observed in the 1997 data with 
an expected background of 1.06 it 0.41 events, and we set an upper limit 

B{Kl -^ TT^e+e-) < 5.1 X IQ-^o at the 90% confidence level. 
13.20.Eb, ll.30.Er, 14.40.A 
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The decay Kl -^ 7r°e"'"e~ is interesting for the study of CP violation and can be used to 
search for new physics. Within the Standard Model there are three contributions |1|] to this 
decay mode. The first is a directly CP violating contribution from electroweak penguin and 
W^ box diagrams; only loops with t quarks, which have amplitudes proportional to rj in 
the Wolfenstein parameterization 0, ultimately contribute. The second contribution is an 
indirectly CP violating contribution [^ due to the Ki component of the Kl- Third, there 
is a CP conserving contribution proceeding through 7r°7*7* intermediate states |Q, which 
can be estimated from measurements f^ of the decay Kl — > ir^'-fy. The total branching 
ratio in the Standard Model is in the range (3 ~ 10) x 10^^^. Branching ratio predictions 
are higher in theories containing exotic particles that contribute to the penguin and box 
amplitudes. In supersymmetric extensions to the Standard Model it is reasonable for 
B{Kl — * vr^e+e") to be as high as 2 x 10~^^, and possibly as high as 10^^*^. The existing 
experimental branching ratio limit is 4.3 x 10~^ at the 90% confidence level. This Letter 
presents an improved limit based on data taken by KTeV in 1997. 

The major components of the KTeV detector, described in detail in reference P,p!^,P|, 



were a magnetic spectrometer for charged particle tracking, an electromagnetic calorimeter, 
and several veto counters to detect particles leaving the fiducial volume. The pure Csl 
electromagnetic calorimeter had an energy resolution for photons of a{E)/E = 0.45% © 
2%/y/E, where E is in' GeV, and a 7r° mass resolution in Kl -^ tt+tt^tt'^ of 1.31MeV/c^. 
The calorimeter was also used to identify electrons by comparing the energy deposit to the 
track momentum as measured by the spectrometer, rejecting 99.5% of all charged pions. 
Additional vr/e separation was provided with eight transition radiation detectors (TRDs) 
located behind the spectrometer P,p!0| . Each TRD consisted of a polypropylene felt radiator 
followed by a two-plane multiwire proportional chamber using an 80/20% Xe/C02 mixture. 
Pulse height readings from the planes were compared with pulse height spectra for pions from 
Kl -^ Ti'^e^u decays that had been identified with kinematic and calorimeter requirements. 
Each pulse was assigned a confidence level for the hypothesis it came from a pion. The 
confidence levels for each of the planes where a hit was associated to a track were combined 



to yield an overall confidence level for the pion hypothesis. For this analysis, cuts giving a 
pion rejection factor of about 35:1 per pion track were used to maximize signal acceptance. 
The trigger required hits in the trigger hodoscopes and in the drift chambers consistent 
with two coincident charged particles passing through the detector. The trigger system 
counted the number of isolated clusters of in-time energy in the calorimeter over ~ 1 GeV 



with a special processor |[Tl]|; at least four such clusters were required. The total energy 
deposited in the calorimeter was required to be over 28 GeV. Events with hadronic showers 
in the calorimeter were vetoed, as were events with activity in the photon veto system. 
Events that passed the hardware trigger requirements were reconstructed online and those 
that passed loose event topology and Csl electron identification cuts were recorded on tape. 

The signature used to look for K^ -^ vr^e'^e" in the detector was two tracks from op- 
positely charged particles with a common vertex that deposited all their energy in the 
calorimeter, plus two extra calorimeter clusters which, if interpreted as photons from the 
charged vertex, had m^^ = m^^o. The decay K^ — * tt^tt^, where vr^ denotes a pion with a 
subsequent vr^ — > e+e^7 (Dalitz) decay, was used to verify the acceptance calculation and 
to measure the number of K^ decays in the data sample. Its signature was similar to the 
signal mode's, but with an additional photon satisfying mee7= ttIt^o. 

A detailed Monte Carlo simulation was used to determine the acceptances for the signal 
and normalization modes, allowing for detector geometry, trigger requirements, reconstruc- 
tion efficiencies, and analysis cuts. It was also used to study background processes. The 
simulation generated kaon decays with the same energy and spatial distributions as the data 
and allowed for particle scattering and inelastic interactions with material in the beamline. 
At the highest kaon energies included in this analysis, a non-negligible Ks component ap- 
peared in decay volume and contributed to the K^ -^ n^-n-^ sample; this was reproduced by 
the simulation and was taken into account in determining the number of K^ decays. 

Final event reconstruction was done after offline calibration of the detector. Recorded 
events were required to satisfy basic quality requirements: the electrons were required to 
not have struck the calorimeter near the beam holes and the track to cluster matching had 



to be unambiguous. The reconstructed vertex was required to lie within the neutral beam 
and well within the vacuum decay region. The reconstructed kaon momentum was required 
the be between 20.3 and 216 GeV/c. Signals used in the trigger were required to lie within 
ranges where the simulation modeled the data well. The tracks were required to be more 
than 1 cm apart at the first drift chamber and have an angle over 2.25 mrad in the lab frame. 
The calorimeter cluster associated with each electron had to have an energy within ±5% 
of the electron's momentum as measured with the spectrometer. Events with four clusters 
(including the two associated to tracks) found by the trigger cluster counter were used to 
search for the signal; events with five such clusters were used to identify Kl — > tt'^tt^,. 

The Ki -^ vr^TT^ sample was selected with requirements on reconstructed mass and total 
squared momentum transverse to the K^ fiight direction (Pj), and had background of 
(0.439 ±0.044)%. Using the acceptance and branching ratio, we determined that there were 
(263.0 ± l.SsTAT ± 15.4sYS ± 9.1br) x 10^ Kl decays between 90 and 160 m from the target 
with Kl momentum between 20 and 220 GeV/c. The systematic uncertainty only includes 
effects not common to both signal and normalization mode, and the third uncertainty is due 
to uncertainties in B{Kl -^ ti^ti^) and i3(7r'^ — > e+e"7) . The largest systematic uncertainty, 
±4.9%, was from the change in acceptance when Monte Carlo events were reweighted to 
match the occupancy in the beam region of the first drift chamber for the normalization 
mode, where tracks are typically close to each other. Additionally, we assigned a ±3.1% 
uncertainty corresponding to reasonable variations in the selection criteria. 

Understanding and suppressing backgrounds is essential for the signal mode search. The 
first background was K^ -^ tt+tt^tt^ where both tt^ showered in the calorimeter; this was 
removed by requiring the reconstructed mass of the event, assuming that the tracks were 
created by pions, to be over 520MeV/c^. 

The second background was K^ -^ ti^tt^ and K^ -^ 7r°7r°7r° with subsequent ir^ -^ e^e~7 
decays. This was reduced by requiring the mass of the two electron system to be over 140 
MeV/c^. There remained some events with two tt^ decays where only one e^ and one e~ 
were reconstructed with a high mass. These events could also have accidentally coincident 



activity. We rejected events with ruee over 370MeV/c^. To ensure that we observed all the 
Kl decay products, we required that the P| be less than 1000 (MeV/c)^. 

The third major background was Kl -^ Ti^e^v with a pion that showered in the calor- 
imeter, accompanied by coincident activity and/or photons radiated from the electron. This 
background limited previous searches |^, but was suppressed here with the TRDs. Tracks 
with TRD data that fit the pion hypothesis with nominally 4% confidence level or less were 
considered electrons; the actual rejection factor, measured using a sample of Kl — * ix^e^p 
decays, was about 35:1 for each pion track. This requirement was applied to both tracks for 
both normalization mode and signal mode; in the normalization mode 94.5 ± 2.2% of events 
passing all other selection criteria passed the TRD requirement. 

Background from Ks — ^ vr^e'^e", assuming that T{Ks — > 7r°e"*'e~) ~ r(i^+ -^ 7r+e^e~), 
is expected to be O(10~^) events. 

The limiting background was the radiative Dalitz decay K^ -^ e"'"e~77 when m^^= 
m.,^0. This background was predominantly internal bremsstrahlung events, but also includes 
Kl — > e^e^7 where an electron radiated a photon while passing through the detector. Both 
components are treated; details are in reference |T^. Figure ^ shows the distribution of 



reconstructed two photon mass, m^^, vs. reconstructed mass of the four particle system, 
rrLee-yy, for data events passing all the selection criteria described above, m^^ was calculated 
assuming that the photons originated at the charged vertex. To the contrary, mee-yy was 
calculated using the 7r° mass and measured photon energies to determine the decay point. 
Use of this "neutral vertex" improved the rriee-yy resolution, permitting a tighter search re- 
quirement on this parameter. The diagonal swath in Fig. |l| is from Kl -^ e~^e~'yy decays, 
where the photons were incorrectly assumed to have come from a vr". 

Three regions are marked out in Fig. [1|. To the left is the region m^^= 135 ± 5 MeV/c^, 
mee77< 465 MeV/c^. Here the otherwise flat distribution of backgrounds from K^ -^ n^e^u 
and Kl — > vr^jTr^, decays with accidental tt^s and from Kl -^ tv^tt^ti'l) decays is peaked due 
to the 7r°s in these backgrounds. This region was not used in estimating the background 
levels. The larger central box in Fig. |l] is the "blind" region, {m^^= 135 ± 5 MeV/c^, 485 < 



^ee77< 510MeV/c^), which was blanked out until the analysis procedure and cuts were 
finalized to avoid human bias. This region was also ignored in the background estimation 
procedure. The smaller central box is the signal region (m^^= 135.20±2.65MeV/c^, mee-yy= 
497.67 ± 5.00MeV/c^) the size of which is roughly a ±2cr range in the signal Monte Carlo. 

To estimate the background in the signal region, we fit the distribution of events to 
{Bq + Bt^ui^^ + Bx'mee'yy + BiH) , wlicre Ti is the distribution in the (mee77,m^^) plane of 
Kl -^ e"^e~77 as determined from simulation, and Bi are the parameters of the fit. The 
estimated K^ -^ e"^e~77 background in the signal region of Fig. || is 36.9 ± 2.0 events. 
The only available method for suppressing this background is to apply phase space fiducial 
cuts [|l^. The optimal cuts were found by minimizing the expected 90% CL. limit on 
B{Kl -^ 7r°e+e^). That expected branching ratio limit was calculated using a hypothetical 
large ensemble of identical experiments in which the only processes contributing to the 
observed number of events were the known backgrounds. The impact of the cuts on the 
signal efficiency, as determined from the simulation, was taken into account. 

The most powerful variables for separating the signal from the background were | cos(^^) | , 
where 6^^ is the angle between the photons and the kaon in the 7r° rest frame, and ^min, the 
minimum angle between any photon and any electron in the kaon rest frame. The pion 
decay angle, which is also the energy imbalance of the two photons, is a good discriminator 
because the comparatively high energy of the photon in a Dalitz decay causes this angle 
to be sharply peaked at | cos(6'^)|=l. In the signal mode, | cos(6'^)| is uniformly distributed 
because the pion is spinless. The minimum electron-photon angle is a good discriminator 
because radiated photons usually appear at a small angle to the electron from which they 
were emitted. In the signal, there is no particular relationship between the photons and the 
electrons and 9^i^ is nearly uniformly distributed. Another widely used variable is |A(0)|, 
the angle between the plane of electrons and the photons. Again, this should be peaked near 
zero for Kl -^ 6^6^77 and uniformly distributed for Kl -^ Ti^e^e^ . After optimization with 
I cos(6'^)| and 6'min, no additional rejection was obtained by cutting on |A(0)|. 

The optimal requirements are |cos(6'7r)|< 0.788 and 6'min> 0.349; with these cuts, the 



expected Kl -^ e+e~77 background was reduced to 0.91 ± 0.26 events with a ~ 24% signal 
loss. The signal acceptance was 3.609 ± 0.017 ± 0.085%, where the first uncertainty is 
from simulation statistics and the second is from the measurement of TRD efficiency in 
Kl -^ tt'^tt^j. The total expected background was 1.06 ± 0.41 events, corresponding (in the 
absence of signal) to an average expected branching ratio limit of 3.3 x 10~^°. 

Several independent checks of the background were made. One important check was to 
verify that the distribution of the phase space variables in the diagonal swath of Fig. |l| 
matched expectations for Kl — > e^e^77 decay. The Kl — ^ e"^e~77 branching ratio mea- 
sured with this sample was consistent with experimental results [0,|ll] and theoretical 



expectations [T^, and the spectra of mee, m^^, and other kinematic variables were well 
simulated. 

Figure ^ shows m^^ vs. mee-y-y for data, including events in the former "blind" region, and 
the 1 and 2 a CL. contours of the Monte Carlo signal simulation. Two events exist in the 



signal region for the data. Allowing for the background [^ and its uncertainty 0, we set 
a 90% CL. limit of 4.85 events. With the above acceptance, dataset size, and background 
level, we set an upper limit B{Kl — ^ vr'^e^e") < 5.1 x 10^^° at the 90% C.L., which is one 
order of magnitude better than the previous result. This limit assumes a uniform three 
body phase space distribution for the signal mode. Alternatively, if the electron pair is the 
product of a single vector meson (as in the penguin and indirect CP violating terms) and 
we model the hadronic vertex with a form factor from Kl — > TT^e^u decay [|l^,|l^, 

/+(g2) = /+(0){1 + [0.0300 ± 0.0016](gVm2)}, 

the 90% CL. limit is < 5.8 x 10"^°. 

If r is the fractional contribution to the branching ratio from direct CP violation, \i]\ < 
4.4A/r at the 90% CL. While not yet as sensitive in probing the CKM parameterization 
as B decays, where recent indirect global analyses ||18| find rj to be below 1, measurements 
in the kaon system are a valuable test of whether the same parameterization is valid for 
both B and K decays. With substantially larger event samples, the use of analyses based 
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on fitting the proper time, Dalitz variables, or decay asymmetries to observed events PJT9[ 
may improve this hmit. 
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FIGURES 
FIG. 1. ra^^ vs. niee'yy for the data after all selection criteria except the phase space fiducial 

cuts designed to suppress Kl -^ 6^6^77. The boxes shown are described in the text; events within 

the central boxes have not been plotted. 
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FIG. 2. As in Fig. |l], on an expanded scale around the signal region, for the data (dots) and 
signal simulation. The contours mark the 1 and 2 a levels for the signal. All selection criteria have 
been applied. 



0.145 



0.14 



> 

O 0.135 



0.13- 



0.125 



• 


• 

• 


• 


• 


• 






• 


• 


( ) 

• 


■( 















0.48 



0.49 



0.5 



nigg^ assuming m^=m^o, in GeV/c 



0.51 
""jtO' 



0.52 
2 



13 



